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Abstract 

A new approach is described to the evaluation of the S-matrix in three-dimensional atom-diatom 
reactive quantum scattering theory. The theory is developed based on natural collision coordinates 
where progress along the reaction coordinate can be viewed as fulfilling the same role as time in a 
time-dependent formulation. By writing the full wavefunction in coupled-channel form it is proved 
that the 3D quantum reactive scattering problem can be treated in the same way as an inelastic 
single-arrangement problem. In particularly, two types of coupled-channel representations, which 
are reduced to two different systems of coupled first order ordinary differential equations describing 
the inelastic scattering, are used. The first system of coupled differential equations is constructed 
on a set of points (grid) of the coordinate reaction curve after solution of many ID Schrodinger 
problems in the directions normal to the reaction coordinate. The second expression for inelastic 
scattering is found using exactly solvable nonstationary ID Schrodinger equation (etalon equation 
method) , which is introduced for describing the localization properties of the full wavefunction along 
the curve of coordinate reaction. In this case we avoid a large amount of computation involved 
in solving the ID Schrodinger problem along the reaction coordinate by using a slightly difficult 
initial conditions for the inelastic scattering equations. In both cases by solving the system of 
coupled first order ordinary differential equations, the full wavefunction and all S-matrix elements 
are obtained simultaneously without further calculations. Our analysis shows that the methods we 
have developed constitute the simplest algorithms for computing the reactive scattering S-matrices. 

Keywords: Quantum Scattering, Natural Collision Coordinate, Etalon Equation Method, 5-matrix Ele- 
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I. INTRODUCTION 



Accurate quantum dynamics approaches have made an immense impact on the theory 
of elementary atom- molecule collisions Recent advances in accurate quantum 

mechanical calculations on simple few atom systems have demonstrated the possibility of 
obtaining detailed information concerning chemical reactions from first principles 4, f|. 
Tremendous progress has occurred in developing and applying both time-independent and 
time-dependent quantum dynamics approaches the last few years. 

At present there are three common approaches used in the quantum reactive scattering 
field: 

1. The hyperspherical coordinate methods [3,0], 



2. The variational methods based on the simultaneous use of mass-scaled Jacobi coordi- 

ll| and 



nates in each of the chemical arrangements involved 



3. Wavepacket methods |12|. Il3j (many references can be found in the review [l^). 

Both time-dependent and time-independent approaches have been successfully applied to 
three- and four-atom reactions. For detailed information at low energy the time-independent 
hyperspherical coordinate approach has proven the most useful {4, 15]. The time-dependent 
wavepacket method has also proved to be very useful, particularly for photodissociation 
and laser-molecule interaction problems. To go beyond four-atom systems, the Multi- 
Configurational Time-Dependent Hartree method has been applied in full-dimensionality 
to calculate thermal rate constants for the H + CH 4 [l6] and O + CH 4 [l7] reactions. 

Despite the successful applications to polyatomic systems mentioned above, general ap- 
plications still appear to be extremely difficult. Quantum mechanics is nonlocal in character 
and numerical efforts to solve the Schrodinger equation increase exponentially with the 
number of degrees of freedom. State resolved accurate quantum dynamical calculations are 
presently limited to at most seven-dimensional problems . Thus, the development of new 
conceptual approaches permitting improved computational algorithms for accurate quantum 
simulations is of great interest. 

Earlier natural collision coordinates NCC allowing an analytical investigation of chemical 



reactivity in collinear collisions were introduced by Marcus 18( and further explored by Light 
[19I . 2(3 ] . Light and coworkers 21, 22] generalized the NCC approach to three-dimensional 



atom-diatom reactive scattering and on this basis they studied H + H2 system and its 



isotopomers in 3D for total angular momentum J = 



Related developments have 



been the formulation and later use of the reaction-path hamiltonian formalism by Miller, 



Handy and Adams 
by Billing 



25 



26 



241 and the reaction path and reaction volume approaches advocated 



27 



28]. 



Due to difficulties in applying the NCC approach, attention turned to other coordinate 
systems. Nevertheless, it seems to us that the primary idea of NCC, namely to simplify 3D 
quantum reactive scattering calculations, can be realized. Particularly, one of the authors 
has shown that for collinear three-body collisions, it is possible to constructed the S-matrix 



exactly, in simple form using the approach 29|, |30(. Seven reactive systems, N 2 + N, I\ 



+ 



3l|- 



O, O2 + O, Li + FH, O2 + N etc have been studied successfully within this framework 

The aim of the current article is to generalize our previous collinear NCC quantum 
scattering approach to a full three dimensional treatment. 



II. QUANTUM MULTI-CHANNEL SCATTERING PROBLEM IN DIFFERENT 
COORDINATE SYSTEMS 



The reactive scattering process in a three-body system can be described by the scheme: 

A + (BC) m , 

(AB) m + C, 

(AC) P + B, 
A + B + C, 



A + (BC) r 



(ABC)* 



A+(BC) m , 

(AB) m + C, 

(AC) P + B, 
A + B + C, 



Sch. 1. 



where A, B, and C are atoms, n and m,p characterize quantum numbers of diatomic states 
corresponding to initial (in) and final (out) scattering channels, (ABC)* denotes an activated 




FIG. 1: Jacobi coordinates for the reaction A+(BC) n — ► (AB) m + C corresponding to the reactant 
(a) and to the product (b) channels. 

complex which may or may not be present as an intermediary species. 



A. The coordinate systems 

We denote the masses of the particles by m^m^ and m>c and the coordinates of the 
particles are denoted by the column vectors ta, i"b and rc describing their positions relative 
to an origin fixed in the laboratory system. The initial 9D problem can be reduced to a 6D 
one by elimination of the center-of-mass coordinates. 

Below we consider the reaction A + (BC) n —> (ABC)* — > (AB) m + C, i.e. there are two 
open arrangements. 

For this reaction reactant and product Jacobi coordinates are illustrated in Fig.l. The 
following equations define the Jacobi coordinates R a and r a (note that index a corresponds 
to (in) or reactant asymptotic channel): 

m B r B + m c r c 
ti a = r A ■ , r a = r c - r B , 

m B + m c 

(R a = K a (R x ,R y ,R z ), T a = T a (r x ,r y ,r z )). (1) 

We now apply the Delves-Smith scale transformation , |32[ to obtain the new scaled 
coordinates qo and qi; 

qo = qoa = A Ra, qi = q la = A^ 1 r a . (2) 



where 



A 



1/2 



m A m b^c 
M 



1/2 



M = tua + rriB + mc- 



(3) 



In these coordinates the Hamiltonian of a three-body system takes the following diagonal 
form: 



1 



#(q;P q ) = — P q + V(g o ,gi,0) ; 



(4) 



where 



q = (qo, qi) = Uk}, fc = 0, ...,5; {%} = (go = |qo|, 9i = |qi|, 92, 93, 94, 95)- 

Note that here and in the following we omit the a channel index for simplicity. In Eq. (jlj) 
fi is the reduced mass defined in Eq. (EJ), P q is the moment vector conjugate to q, in the 
body-fixed system and 9 is the angle between vectors qo and qi (see FIG. 1 and FIG. 2). 
An alternative coordinate system (go, 9i, 92) is illustrated in FIG. 2. Either of the coordinate 
systems (go, gi, Q2) and (go, gi, 6) can be used to describe the motion of body-fixed system in 
the plane of the vectors (qo, qi) and they will be called intrinsic coordinates. The interaction 
potential between all atoms depends only on the intrinsic coordinates. 

The coordinates (q^, q±, q$) are taken to be the Euler angles fl(= , (2, G), which orient 
the three-body system in the space-fixed frame (see for example refs. 13[, |33f] and 34l|). 

The optimal coordinate systems for reactants and products are different jo), [3^]. This 
fact creates certain mathematical and computational complexities in the investigation of 
the reactive scattering problem. One way to proceed is to turn to a special type of curvi- 
linear coordinates which continuously leads from reactant channel to product channel. The 
possibility of introducing the coordinate system which is simultaneously suitable for descrip- 



tion of (in) and (out) asymptotic states was first time discussed in the work [36|. At last 
corrected version for collinear collision was presented in the work 18J and named natural 
collision coordinate (NCC) system. 

In order to define the NCC appropriately, we connect the (in) and (out) asymptotic 
channels by a smooth curve Sj/. The curve depends on two parameters, one of which can be 
chosen arbitrary and the other is system dependent. As progress along the curve represent 



the passage from reactant to product during this it will be called reaction coordinate 



19]. 



B 



FIG. 2: The three-body system in mass-scaled Jacobi coordinates. The coordinate go shows the 
distance between particle A and the center of mass (BC) pair O, the distance between particles B 
and C is q\ and q2 correspondingly between A and C, 8 is a Jacobi scattering angle. The point 
O is a center of mass diatom (AB) and 9 is correspondingly Jacobi scattering angle which is used 
in (out) channel. 

The curve 3^, which connects (in) and (out) asymptotic channels can be determined in 
the plane (go, gi, # = 0), by simple formula: 

<fo = c a _ + K<?1 C - <Uq) + <lt q , <hq < Ql < +OC> (5) 

1l ~ Qeq 

where a and b are constants. 

Note that the collinear collision configuration A+(BC) describes on the plane (qo,qi,9 = 
0) while the collinear collision A+(CB) describes on the plane (go, Qi, = tx). The motivation 
of choice of curve for later definition of 3D NCC system is given in the Appendix A. 
In Eq. §5§ q~ q and g+ are mass-scaled equilibrium bond lengths of molecules in the (in) 
and (out) channels correspondingly. Note that a is an arbitrary constant, which is usually 
chosen to make the curve pass close to the saddle point of the reaction. The constant b 
depends on the mass of particles and is given by expression b = cotcp = \mj<jnc / m BM'\ 1 ^ 2 . 
The superscript c over go and gi underlines the fact that the point (gg, qf) lies on the curve. 

The limit gj — > q~ q , q$ — > oo corresponds to the reactant region or (in) channel 
while the limit ql — > oo, q^ — >• oo on the curve 9^ corresponds to the product region 
or (out) channel, where gg = bq\ + g+ It is easy to show that the mass-scaled dis- 



tance q 2 between A and B particles connected with the curve 3^ is described by ex- 
pression q% = 92 L — \J (Qo) 2 ~ 2^9o9i cos ^ + ^ 2 (9i) 2 5 which in the (in) channel is equal 
{q.2 ~~ 9o) Igc^q- = —bq^qCosO. In other words it is denote that in the (in) channel curve 9*^ 
is coming to the equilibrium distance q~ q of (-BC) diatom. In the (out) channel the curve 
approaches the equilibrium distance g+ of diatom (AB) (see Appendix B). Note that 
the Jacobi scattering angle 9 which defined in the reactant (or (in)) coordinates system in 
the product (or (out)) channel obviously must be limited to inasmuch as the two sides 
BO = bqi and OA = go of ABO triangle are limited to infinity while the third side is lim- 
ited to finite quantity AB = q 2 — > q~L. Remaind that it doesn't mean that the rotation of 
new-formed diatom (AB) stops. The rotation of diatom (AB) is described by angle 9' (see 
FIG. 2) which in the reactant region is limited to but beginning from strong interaction 
region up to (out) asymptotic channel obviously can have any value. 

The coordinate u describes the translational motion of three-body system between re- 
actant and product channels and is changed along the curve measured from an initial 
point Uq. It particularly can be determined by equation: 

9l q.eq 

The signed distance from the curve 3^/ in the plane (qo,qi,0 = 0) is given by v. We can 
organize the one-to-one mapping between coordinate systems (go, Qi, &) and (u, v, 9) in only 
some subspace of internal 3D configuration space. It should be emphasized that above 
mentioned subspace must include a part of 3D space, in which the probability current of 
reactive scattering process is localized. This condition can be fulfilled if the coordinate 
reaction curve is correctly determined for a configuration of reactive collision with the 
minimal energy surface (see FIG. 3). It is obvious that in this case the condition of one-to- 
one mapping between coordinate systems for any configuration of reactive collision can be 
automatically satisfied. 

After satisfying aforementioned condition it is easy to write transformations between set 
of coordinates (<Z0)9i) and (u,v) jsj]]: 

qo(u,v) = q c Q (u) -vwa.<j){u), 

qi(u,v) = q{(u) + v cos0(u), (7) 
where the angle 4>(u) is determined from the requirement that the coordinate system (u, v) 




FIG. 3: a) Intrinsic 3D Jacobi and natural collision coordinate (NCC) systems. The angle ip is 
defined in the text, b) The reaction path (cycled curve) is passing through the minimums of potential 
energy while the reaction coordinate 9^ can be an arbitrary smooth curve connecting (in) and (out) 
asymptotic channels. The lower shaded area is self-crossing region for NCC system associated to 
the reaction path curve while the upper one self-crossing region for NCC system is associated to 
the curve Gjj. 



should be orthogonal (see Appendix C): 

= cot <ft(u), lim cot0(w) = cot tp. (8) 



(u,v=0) 



In the NCC system the motions of the body system are locally factorized into translational 
(w, infinite extension), vibrational (y, finite) and intrinsic rotational (9, finite) motions. 
Note that the coordinate u is perpendicular to the curve v and independent from the angle 9. 
These properties ensure that the full wavefunction can be conditionally (locally) factorized. 

So, we defined NCC system for investigation of 3D reactive quantum scattering which for 
using has only one limitation related with the region of self-crossing of coordinate lines. In 
other words it is necessary to choose the reaction coordinate in such a way in order to 
the probability current of quantum reactive scattering process in the region of self-crossing 
of coordinate line was absent. 



III. QUANTUM REACTIVE SCATTERING IN THE NCC SYSTEM 



A. Equation of motion of the three-body system 

The overall stationary wavefunction of three-body system can be written as: 

* JM (q , qi ,9,n) = —J2^( ( lo,qi,0)D J KM (n), (9) 



qoQi K 



where J is the total angular momentum, M and K are its space-fixed and body-fixed z- 
components respectively. The summation over K ranges from — J to +J in unit steps and 



D J KM is the Wigner D-matrix [33, |34| . 



After separation of the external rotations, the action of the hamiltonian on the wavefunc- 
tion is given by the formula 0| : 



A j h 2 { d 2 d 2 1 —j 

V mi J I sm 



h 2 fq 2 + q 2 \ f 1 d . d K 2 \ j 



2 /A mi J {sm9 89 89 sin 2 

»-2 

+ V(R, r, 9)$ J K (q , q u 9) + —^ (j{J + 1) - 2if 2 ) $ J K (q , q u 9) 
h 2 [ 8 1 — j 

h 2 [ 8 1 — j 

+ 27^f CjK \M + (A ' " l) cote j $ (^-i)( go ' 91 ' e) 



E<f> J K (q Q , qi ,9), (10) 



where 



[J(J+l)-iT(iC±l)] 1/2 . (11) 
After the coordinate transformation (qo,qi,9) — > (u,v,9) in equation ( flOl) we find: 



H$ K (u,v,9) 



2jj,r]\du r] 8u dv dv J K 
K 2 Uj 



/? 2 



2/^ g 2 (w, u) 



x 



1 g • 
sin^^ Sm <9# sin 2 fl 



$ k (w,M) + ^(«,M)$kKM) 



+E JK $ J K (u, u, 0) + — <j H+ K ® J {K+1) {u, v, 9) + HjK&lx.q {u 
E<5 J K (u,v,9), 



;i2) 



where 



E JK (u,v) 



h 



2(iql 



— (j(J+l)-2K 2 ), H 
u,v)\ / 



C ^ K (K± l)cot0 



q$(u,v) { 89 



q 2 (u,v) 



Qo(u,v)qf(u,v) 



V(R, r, 9) = U(u, v, 9), $ J K (q , q u 9) = & K (u, v, 9). (13) 



q%(u,v) + q 2 {u,v)' 

In expression ( jl2l) the Lame coefficient rj(u,v) has a form (Appendix C): 

(is 

77(11,0) = [1 + , 



(14) 



where K (u) is described curvature of the reaction coordinate 9^ in the continues point w, 



correspondingly s the length along the 9^-: 

2a 



iT(w) 



6 - a/ (q£ - g e(? ) 



-3/2 



and 



1 + 



1/2 



b-a/fa-q^) f s6 + a/(5i-9j 



(15) 



(16) 



Finally we transform Eq. ( 1T21 by setting <& K (u,v,9) = rj 1 / 2 $^(m, w, 6 1 ), whereby we 
obtain: 



did d 2 



du t] 2 du dv 2 



q 2 (u, v) 



1 d . n d K 2 



sin9d9 89 sin 2 # 



+ ^[E-U(u,v,9)- 



E JK (u,v) -U eff (u,v)\ }® J K (u,v,< 



{H+ K $ J {K+1) (u, v, 9) + Hj K * J {K _ x) (u, v,0)j = 0, 



(17) 



where the effective potential U e ff(u,v) is defined by 



1 d 2 r] 5 / dr] 



(18) 



4?7 2 \ dv J 2r/ 3 du 2 4?7 4 \ du { 

Schrodinger equation ( 1TT1) is in a form which is suitable for further investigation by the 
coupled channels method. This will be done in section III. C after the 5-matrix has been 
defined in Section III. B. 



B. Definition of 5-matrix in term of reaction coordinate 



In the standard scattering theory the main problem is the construction of the scattering 
S'-matrix, holding the transition amplitudes. Let us discuss the exact representation for the 
S'-matrix in terms of overlap between stationary wavefunctions (see for example {ssl, 3o|). 
In the body-fixed NCC system we can write the following general formal expression: 

n'j'K' 

where $„^( u a , v a , 9 a ) and ^ n ijix\ u Pi v Pi @p) are the total stationary wavefunctions in cor- 
respondingly NCC systems which are evaluating from some clear (in) and (out) asymptotic 
states, n and j are vibration and rotation quantum numbers, S^,j, K ^ n j K are an complex 
matrix elements. 

For the overall wavefunction $^k(u,v,9) we may enforce the following asymptotic be- 
haviors and boundary conditions: 

n'j 'K' 

lim $S J (u, u, 9) = lim ^-^(u, u, 0) = 0. (20) 

where (m, #) and $i°"/^/(w/3, U/j, ^) are an asymptotic wavefunctions correspondingly 
in the (in) and (out) channels, S J n j K ^ n ,j, K i are a reactive S -matrix elements. 

Similar conditions can be also written for wavefunction &^jx(u, v, 9). For later calcula- 
tions it is important to define the behavior of the wavefunction (u,v,9) in the (out) 
asymptotic channel: 

^( U p,Vp,9p) u T+oo ^ J (u Pt vp,0 p ). (21) 

Note that in fl20l and fl2Tj) the functions $^™^ J (w, t>, 6>) and &njK J ( u Pi Vp, 9p) are described 
by expressions: 



= -^=exp(-ip- jU )U^(v)Q jK (9), 
= exp (ip+ up) (v p ) e jK (9,3), 



^\up,Vp,9p)^^exp(ip+ U p)U^(vp)e jK (0p), vl 3 = \ sJME - e^), (22) 



where e is the vibration-rotational energy of the initial state of the diatomic, is 
vibration- rotational energy of the final state, Qjk{0) is a normalized associated Legendre 



polynomial 



34|, I&?\(v) is the vibrational wavefunction of the initial diatomic, which sat- 



isfies the following equation: 

h 2 d 2 



2/i dv 2 



U {in \v) + 



fr 2 jU + 1; 

2fiv 2 



e~ Tl (in) (v) 



(23) 



where U^ n \v) is bounding potential of diatom. Similar expression is valid for IL^f\(v) after 
replacing U^ n \v) — > U^ out \v) and e~,-^ — > e t,(jy 

The matching conditions between a and f3 full internal wavefunctions can be write in the 
following form (see for example fill, iiol]): 

(24) 



K 



K 



The relation between the rotation 
know may be represent by formula 



unctions in different body-fixed frames in (1241) as well 



11 



40 



4l|: 



^Mi?(Cl/3; C2/3, C3/3) — ^ d J MK ,( r d)D'^ K ,(Ci a , ( 2 a, (3a] 



(25) 



K' 



Note that d is the angle between vectors qo a and qo/j, and may be described as a function 
of the internal coordinates using the relation: 



qo Q qo/3 = % a qop cos 



(26) 



where $-angle between a and (5 body-fixed reaction coordinates systems. 

Using the expressions (1241) and (1251) it is easy to find the following equation of connection 
(Appendix D): 



<f>^(u a ,v a ,9 a ) = J^C^C^'^'^ 

V ' P 



(27) 



A 



Now taking into account that the overall wavefunction ^^)^(up, vp, 8p) at the limit u 



+oo going to clear (put) asymptotic state ( 1211) from equation ( Tl9|) we can fined the following 
expression for scattering matrix (see Appendix D): 



(...) 



u0 



1/2 dv / sin0d0. 



(28) 



Putting the expression for $^(« a , v a , 6 a ) from (j2"T|) in the ff28|) we can find: 

Sn'j'K' ^njK(E) = 



^ lim 

rt . u—>+oq \ J J 



v6 



(29) 



In the expression fl29|) S n ijiK'^ n jK maybe interpreted as a general form for the reactive S- 
matrix elements. Evidently from definition of the S-matrix elements ( |29l) the integration over 
coordinate u is absent, however it is easy to understand that in the (out) asymptotic region 
the full phase is cancelled out and thus S n ijiK'^njK -matrix elements become independent 
from coordinate u. 

Now the differential reactive state-to-state cross section may be simply constructed with 

(E) (see particularly works Q, Q, 42]): 



the help of matrix elements S n ,j, 



n'j ' K <— njK 



(30) 



KK' ™ 3 J 

where p n j is the incoming wave vector. 

The total integral reactive cross section for reaction from a particular initial state to all 
possible final states is then given by a summation over all total angular momenta which can 
contribute to the reaction: 



a 



tot 

all*— nj 



V 



£(2J+l)P/(£), 



7T 

n i j 



where 



P r J (E) 



-njK 



(E) 



(31) 



(32) 



n'j' KK' 

Recall that P^(E) total reaction probability for a particular value of the total angular 
momentum J. 

So, we have obtained a new representation (129]) for the scattering S'-matrix, where the 
coordinate u, which varies along the reaction curve 3^/, plays a role much as usually time 
does in standard quantum scattering theory. The reaction coordinate u could thus be viewed 
as an intrinsic time of the scattering. 



C. Coupled-channel expression for full wavefunction and 5-matrix elements 



Remembering that along u the system is in translational motion, while along v and 6 the 
motion is localized we employ the time-independent coupled-channel (CC) approach. 



It is convenient to write down the intrinsic full wavefunction in the following form (see 



for example 



(33) 



"j 



and necessitate the wavefunction fl33l) to satisfy to initial conditions fl20l) . the symbol [g] 
describes a set of initial quantum numbers. Note that it is standard way to obtain the close 
system of coupled equations for scattering (or translational) functions G 



e'[e] 



u . 



The vibrational part of the wavefunction , E n ^(u,v) forms an orthonormal basis in the 
variable v for each fixed value of u and satisfies the equation: 



h 2 d 2 - ( . h 2 j(j + 1) 

- + U(u,v) + — ^ 



2/x dv 5 



2/jlv'' 



l n (j)(v] u) = e nij) (u) E n{j) (v; u), 



(34) 



where 



U(u,v) = U(u,v,9)\ -U eff (u,v). 



(35) 



Recall, that U(u, v, 9) \ _ Q corresponds to the potential energy of collinear collision. 

In some situations it is useful to approximate U (u, v) such that Eq. (134)) can be analyti- 
cally solved (see section E). Note that the (}3~4l) in the limit u — > — oo transforms exactly into 
the asymptotic equation (1231) . Correspondingly in the limit of u — » +oo Eq. (|34p describing 
bound state of the (out) asymptotic channel. Note that for both cases u — > =Foo the effective 
potential energy U e ff(u,v) limits to zero. 

For subsequent analytical manipulations we give two important expressions for Legendre 



polynomials 



33 



1 3 d 
sin 6- 



K 2 



sm9d6 d6 sin 2 6 



e jK (9) = - 3 (j + i)e jK (9), 



(36) 



and 



Q j{ K±i } (6) = cf K e jK (6), c% = [j(j + 1) - K(K ± 1)] 1/2 . (37) 



T^ + (K±l) C ot6 



Next substitute the wavefunction expression (1331) into Eq. (fl7|) taking into account (I34|) - 
( |37l) and multiply by E n '^>)(v;u) and Qj>K'(0)- Thereafter integrate over the angle 9 and 



coordinate v to find the following equation: 

d 2 / <9 ldrj\ d t I d 2 2d?] d 



n du 2 \du rj du ) n , n du \ <9u 2 77 <9u <9u / n ,_ ~^ h 2 ^ 



2/i 



£ - E JK r(u,v) 



- e n(j)(u) + U(U,V) + 



2uv 2 



+ <^2 / ,Jj'J 
C[q I n'n 



= 0, (38) 



^Jj K'^fij (ftT'+l) + ^Jj K'^nj'iK'-l) [e]^ 



where Ct^ x = Cjk^k ( see expressions (TTTT) and ([3] 

The equation (138B may be presented in another form: 



d 2 2 l—--?^\ 2/x 
1 n <iw 2 \ du T) du I , n du \ du 2 r\ dudu / ,_ ft 2 



- e n(i)( u ) + U(U,V) + 



2fw 2 



E - e jk(u,v) 

^K'K 



(+)/ 

nj K [g] 



u) = 0, 



(39) 



where the summation over repeating index n' and j ' are implied and we use the following 
notation for matrix elements: 



(f( U ))nn> 



U«,(u,v) 



E n(j)(v] u)f(u, v)E* n , {j) (v] u)dv, 

o 

Q jK (0) U{u, v, 6) & fK (6) sm6d6. 



(40) 



Thus we obtain a system of N coupled second-order differential equations (!38|) or (1391) . 
which can be rewritten in a form of 2N coupled first-order ordinary differential equations 
system: 



du 



I 
A B 



(41) 



where 



\F(u) 



i<N 



G { +)?Au) 

ri] K[q\\ > 



i>N 



dG[%, ,(u) 

n 3 K [gj v y 

du 



. 1 . 



(42) 



Moreover in Eq. (|4"TI) the following denotations are made: 
d 2 2 dr) d \ 1\± 



A n >j'K' \ nj K{ u 



du 2 r]dudu/ n , n ' h 2 \^ 



V 



E - E jR (u,v) - e n{J) (u) + U(u,v) 



5j ,j 5 K 'k 



(/q / nrr 



Sk'+1K C jJ(K-1) + 5 K'-1K C Jj(K+l) 



Id 1 dr\ \ 



(43) 



Now we turn to the derivation of S-matrix elements using the representation for the full 
wavefunction fl33|) . 

As it is known the full wavefunction of three-body system is determined by the set of three 
quantum numbers (n,j,K). If we want to describe the full wavefunction of body system 
with the help of representation (13"3"|) it is necessary to demand the asymptotic condition: 

1 



lim Yd?" 

M -no n'j 'K' 



U) 



exp (-ip nj u) b nn <bjjib K K< 



(44) 



to be fulfilled. It means that in the limit u — > — oo, the wavefunction (133]) transforms to the 



asymptotic wavefunction $ {u, v, 6) (see formula (j20j) ). 



The 5-matrix elements are obtained after substituting (133]) into the expression (1301 ) : 



Sn'j'K'^njKiE) = \r^=- ^hm [(?] («) WW («) Aj^/ 



(45) 



where 



Wn'n(«) 



; R(j)( u ; w ) n SSV v ' 



A-jK'^jK' = (e 3K ,(9)Q jK ,(9)) = 5 3r (46) 



The expression for ^-matrix elements (I45p can be simplified, if we take as basis the functions 
E n (j)(v; u), which in the limit u — > +oo coincide with the orthonormal basic wavefunctions 
Tr^)(v). In this case we get the simplification \im u ^ +00 W n in(u) = 5 n / s , and the following 
expression holds for S'-matrix elements: 



Sn'j'K'^njK(E) 



P n'j' r (+)J( W . ,„<, 



P 



q = (njK). 



(47) 



So, we have now shown that the initial quantum scattering problem can be exactly 
reduced to the standard set of N coupled second order differential equations fl38l) - fl39l) in 



a single variable u (where N is the number of channels or coupled equations). For the 
solution of this system of equation it is useful to represent it to the form (|41|) . The matrix 
equation (141 j) has 2N linearly independent solutions. There are several standard ways in 



which this set of equations may be solved 



45 



46]. The 2N independent solutions are 



then combined to give N solutions which obey the asymptotic conditions as set out in Eqs. 
(JHJ). The process of combining the independent solutions to satisfy the boundary conditions 
automatically yields the S-matrix elements. 



D. Etalon equation method, full wavefunction and 5-matrix elements 

During the development of the algorithm for numerical simulation based on aforemen- 
tioned theoretical approach it is necessary to recover the scattering matrix equation (I41j) - (]43p 
on the grid. The grid is formed along the reaction coordinate curve, which connects (in) 
and (out) scattering channels. The problem is that at each point of the grid one has to 
solve Mi = iii x ji number of one-dimensional quantum problem and by further numerical 
integration to find the form of the matrix equation (f4"T!) on the grid. This procedure re- 
quires huge amount of computer resource and undergoes the accumulation of computation 
errors. Besides matrix equation (T4T!) in this case will be given as a numerical array, which is 
extremely inconvenient for numerical simulation with changing integration scale. In order 
to overcome this difficulty one can use the exactly solvable model for vibration states in 
the total wavefunction coupled-channel expression. This allows getting analytical form for 
one-dimensional matrix equation of reactive scattering. 

So, for full wavefunction can be written another coupled-channel representation: 



= E47km,]N S ?(J)(«^)%^W» < z7 < oo, (48) 



v .1 

where s[ s ,"uu, v) is the wavefunction of singular nonstationary (on reaction coordinate w, or 



later intrinsic time r) quantum oscillator and satisfying the following equation 47J (later 
will be named the etalon equation): 



Q£ (so) 

ih 



S M E^ so \t;v) = E^ so \u,v), 



d 2 f ^u 2 (t) 2 g_ 

dr [ ~2udv 2 V 2 v 2 , 

h 2 

r = u/vq, g = —j(j + 1), < v < oo, — oo < u,r < +oo. (49) 



Note that vq = a^Ui (where a is a Bhor radius, for explanation u,i see (153]) ) give the 



characteristic speed by coordinate of loca 
The Eq. 



can be solved exactly [48j] for values v > 0: 



ized motion v. 



JL 

fur' 



a+l 



Tin 



1/2 



a 



T(n + a+ 1) 
da . G?7 

dr' dr 



,0+1/2 



2V- ff7 • (5 °» 
It is significant that the solution (]50|) may be analytically continued to the infinite numerical 
axis —00 < v < +00 (see particularly (|47j])). In (150]) the function 7(7-) = J u \a(r)\~ 2 dr is 
the argument of complex function a which is a solution of classical oscillator problem: 



a + us r cr 



0. 



(51) 



For further investigations it is useful to have an analytically solvable model for function 
cr(r). In particularly the Eq. (|5"T1) may be integrated exactly if the frequency is described 
by following model form (see Appendix E) 491 ] : 

uj 2 (t) = A + A 1 tanh(Ar), (52) 

where Aq, A\ and A are some adjusting constants, which can be chosen to be most suitable 
for numerical simulations. As it is evident from formula (152]) . lim UT _ 
corresponding asymptotic solutions of Eq. (]5"T]) are: 



>=FOG 



ufj and 



<y t 



o- r 



r— >+oo 



- 1/2 K^l 
-OO c 



10 



-1/2 



where Uij are coefficients at second terms of expansion series in corresponding asymptotic 
bound states potential energies (see Eq. (|23l ). The coefficients Ci(2) are some complex 
numbers which are found by the solution of Eq. f]5"T]) in the limit r — > +00. The constraint 



A T A, 



(53) 



\Ci\ 



I Co 



1 is hashed at the numbers Ci and C2, which is followed by commutation 



correspondence. 

Using this fact it is easy to calculate bound-state energy in asymptotic channels e v (j) = 
fvjjij{2v + a + l) and the asymptotic form of the etalon wavefunction (150]) . In particular, in 
the (in) asymptotic channel: 



h 



a+l 



T(n + 1) 
T(n + a + 1) 



1/2 



> a+1/2 exp(-2ivUiT 



v 

2h 



(54) 



In the Eq. (1541) L a v {hr x u'jv 2 ) is Lager polynomial and the wavefunctions (r; v) form 
an orthonormal basis: J °° 5[!°|(t; t?)Hw°l*(r; t> )c/t> = 

When g = (a = 1/2) the wavefunctions (l54"j) describe the nonstationary oscillator 
with wall at the beginning of coordinate 50j. It is well know that the wavefunction of this 
system coincides with odd functions of harmonic oscillator in view of identity H2 V +\(x) = 
(—iy2 2u+l v\xL l J 2 (x 2 ) (see for example 51])- Nevertheless the Eq. (H9l) have another type 
solution too, which coincides with even wavefunctions of harmonic oscillator 52]. 

So, in case of g = the solution of etalon equation (j39]) is wavefunction of q uantum 



harmonic oscillator depending on time (in this case intrinsic time r) frequency [43, [53] : 

v/2 



■=:( s °) t 



T-V) 



a 
27 



(z/lav^r^exp^ 2 )^^ 



(55) 



Now comparing representations ( J34l) and ( 1481) is easy to find the following relation between 



reactive scattering functions G 

G 



(+)J 

n]K> [g] 

(+) J I \ 
Ut 



[u) and G { °? ] K i [e] (u) 



u)G 



(so)J 



v 

where matching coefficients C m Q^(u) is defined by formula: 



u 



(56) 



Cpfi(j ) [U 

Using (l5Tj) in the limit u — > — 

■y(so)J 



(57) 



oo for definition of asymptotic conditions for the function 



G_^°J,,_t(u) the system of linear equations can be obtained: 



^i>0(j)\ a ) ^ujK'[g] 



(u) = 0, 



E- c 



i(so) J 



(58) 



Y (S0)J 



where m defines the number of vibrational states in the (in) channel in the same way as 



ol u 



) = (2tt) 1/2 exp(-ip nj u) and C. 



.+00 ,-,(-) ( 
-oo 



uv 1 ;v)E* n{]) (u,v)dv. 



In the system (l58~I) the number of equations is finite while the number of unknown quanti- 
ties is infinite. For exact definition of problem (1581) we have to demand following conditions: 



G^niu) 
vj K' [g] V J 



0, v > m. 



(59) 



Now substituting (|48|) into (fTTj) and after simple analytical calculation we can get a new 
equation for reactive scattering: 



djij o K i K 



d 2 Id ldr}\ d 



d 2 2 ( ldr] I fv } — 
du 2 \du 7]du/ u , D du \du 2 \rj du h~ J du / v , t 



E - E jR (u, v) - U eff (u, v) + ^ V + ^ (j + l > 



, % I v'v 



$K'+1K ^jjK-l ^K'-1K C 



JjK+l 



(60) 



In this case the summation over repeating index v is implied and we use the following 
notation for matrix elements: 

/+oo 
Et u) (uv^;v)f(u,v)E s 1/ f {j) (uvo 1 ;v)dv. (61) 
-oo 

Remind that the system of second order differential Eq. (I60p may be written in the first 
order form of type (l4Tj) -(T43l). 

Finally, taking into account (H^l)-()46l) (HHl) and (1571) one can find the new expression for 
^-matrix elements: 



S n 'j'K> ^njK(E) = J lim ^ G^, [g] (ll)C m Q) (u) W n n> (u) X jlK , _q K , 

Pnj u ^ 



IPn> 



P 



un'(j 



(62) 



"J 



So, we found a new analytical expression for the S -matrix elements of reactive scatter- 
ing (l62l) with the help of exactly solvable etalon equation fj49|) and we will call it etalon 
equation method. Obviously this method is favorable for numerical simulation. Particularly 
computation of S-matrix elements (|62j) in this case is relatively simple too. 



IV. CONCLUSION 



The introduction of natural collision coordinates by Marcus 18| was intended to simplify 
quantum reactive scattering calculations. In spite of significant efforts [19|, |20|, |2l| the 
application of this method to 3-D reactive scattering has encountered considerable difficulties, 



and as a consequence the investigations in this direction ceased in 1988 (see for example 
report Q]). In particular, the following two problems can be indicated: 

1. For an atom-diatom arrangement reaction A + BC in 3-D, at every conserved total 
angular momentum J, the vibrational coordinate v (see FIG. 3) becomes a 2D (ro- 
vibrational) surface in the 3D intrinsic space which specifies the size and the shape 
of three-atom triangle. It is significant that this surface has extremely complicated 
the metric properties and can be investigated with the help of difficult numerical 



calculations 



22)- 



2. The third arrangement reaction AC+B (see schema 1) must also be included in the 3D 
computation schema. However, this leads to problems in the understanding of trans- 
lational reaction coordinate idea. Moreover, technical difficulties arise by matching 
surfaces between the three arrangement channels. 

Nevertheless the mentioned difficulties either can be overcame in the intrinsic 3D space or 
they do not play an essential role in the framework of this consideration. 

Hence, if we define three different curves of the reaction coordinate analogously to 
§5§ by the way they connect different asymptotic subspaces ( different reactant and product 
channels) one can describe all the scattering channels (see schema 1). From the other side, 
by defining this curve and correspondingly NCC system carefully on the case of collinear 
collision all possible configuration of reactive scattering may be described. 



In recent articles 



29 



30 



3l|, one of the authors together with colleagues have analyzed 
in detail the difficulties arising in the formulation of quantum scattering theory in NCC 
system. 

In particular, it was shown for the case of a collinear scattering process in NCC system 
that a procedure similar to that described above for the 3D case could effectively yield the 
exact S-matrix. 

In this article the quantum reactive scattering theory in curvilinear reaction coordinates 
has been generalized to the 3D case. 

We have started with Schrodinger equation of thee-body system in Jacobi coordinates 
(TTU]) . Using the coordinate transformation given in Eqs. ©-flED this equation has been 
transformed into the NCC system. The reaction coordinate u (or r) may be considered as a 
chronological parameter in the theory (i.e. intrinsic time), and plays a role analogous to that 



of time in standard time-dependent scattering theory. Representing the full wavefunction 
of three-body system by standard coupled-channel form fl33|) the initial 3D multi-channel 
quantum scattering problem (TTTj) with correspondingly initial and border conditions (120j) 
may be reduced to an inelastic single- arrangement problem Eq. fl5Sl) - fl591) or fHTjl - fHS]) . As 
a result of this reducing after the solution of system of Eq. (HTi) - (H5l) on a set of points u 
of curve of reaction coordinate 0^/ (grid) the full wavefunction and all transition S-matrix 
elements are found simultaneously (see expression for S-matrix elements fjTT|) ). 

Another direction of investigation of Schrodinger Eq. (fl7|) is the representation of solu- 
tion for full wavefunction of three-body system in the coupled-channel form ( 1481) combined 
with the exact solvable etalon equation method (see Eqs. ( I49I) - (15T)|) ). Recall that etalon 
equation describes the localization properties of full wavefunction along the curve of reac- 
tion coordinate 3^/. This method reduces the 3D quantum multichannel scattering problem 
to inelastic single-arrangement problem too. Still, the initial conditions for the system of 
differential equations ( |60l) in this case are a little bit complicated and can be found by the 
solution of linear algebraic equations system (l58l)-(|59l). 

Note that this is an important theoretical result, which seems to very useful for numerical 
calculations. Remind that in traditional approached the ID scattering problem Eq. fl34|) - 
(|35|) is constructed after very grate volume of grid computations of ID Schrodinger problem 
along the reaction coordinate. 

The main theoretical advantage is that the wavefunction and all the body-fixed S-matrix 
elements are determined from the solution of standard coupled differential equations in only 
one variable u (the scattering coordinate) simultaneously. The body-fixed S-matrix is used 
to determine all the differential and integral state-to-state reactive scattering cross sections 
of the system (see formulas fl30l) and fl3TT) ). 

It is obvious that on the basis of developed S'-matrix representations the maximally- 
possible effective parallel algorithms for direct numerical simulation of 3D reactive quantum 
scattering problem may be elaborated. Particularly numerical solution of both systems of 
inelastic scattering equations ( 1391) or (T411-(l43l) and ( 160|) can be realized with the help of the 
i2-matrix propagation method simultaneously yields the full wavefunction and all S'-matrix 
elements without further calculations. 

Finally, we would like to note that the usage of NCC system in quantum multichannel 
scattering theory will permit to carry out similar type reduction for any amount of atoms. 
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FIG. 4: T/ie curve of reaction coordinates for two different exchange reactions on a collinear colli- 
sions plane (R,r). It is obvious that in this plane there are not interferences between two exchange 
reaction channels. Recall that angles ip\ and (f2 are given by formulas cot(fi = [mAmc / m BM~\ 1 ^ 2 

1 /2 

and cotif2 = [mAmB/mcM^ 

This will essentially simplify the calculation. 
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APPENDIX A: CHOICE OF 3D NCC SYSTEM 

Two collinear collision configurations between free particle A and bound pair (BC) 
schematically may be represented on the plane of Jacobi coordinates (R, r) in the following 
case (see FIG. 4): The smooth curves (coordinate reactions) 3^/ and 3^/ in figure con- 
nect correspondingly reagent and product asymptotic subspaces in the collinear collisions 






FIG. 5: Energy surfaces of reactive system H + H 2 in the plane of Jacobi coordinates (R, r) for 
different angles 6 = 0°, 30°, 60°, 90°. As it is evidently in the plane a) two reaction channels are 
opened while on the planes b), c), d) only excitation channels are opened. Note that it is a universal 
property which follows from kinematic geometry of scattering and doesn't depend on what A, B and 
C particles are interacted. 

A + (BC) and A + (CB). It is obvious that in this case the different reaction channels are 
isolated and there is no interference between them. The illustration of this fact is shown 
on the example of ah initio PES, termed LSTH, type potential energy surface for reactive 



system H + H 2 (see FIG. 5a) [56|, |57|. When the Jacobi angle is fixed 6 £ (0, re), going along 



the curve 9^/ or 0^ we can't come to mentioned products subspaces. In other words the 
exchange reaction channels on the plane (R, r) in this case are closed and this is an universal 
kinematic property of scattering and doesn't depend from particles sort (see FIG. 5b-5d). 

However, in this case the reaction goes out of mentioned plane in 3D space (see FIG. 
6a-6d and FIG. 7a- 7d). As kinematic analysis shows, the 3-D region, where the probability of 
reactive scattering processes is concentrated, is essential around the collinear collision plane 
(R,r). Moreover, this region is compressed to the plan (R,r) as system is evolving along 
coordinate reaction curve Q 1 ^ (or Sfr'j/) to the asymptotic subspace of products. Recall that 
all topological properties of reaction surface, which are in Jacobi coordinates system, are 
conserved in 3D NCC system too. 

So, the 3D NCC system, which is connected with the reaction coordinate curves 9^/ 
or Ss'if, describes uniquely 3D part of configuration space, where multichannel scattering 
process is concentrated. 




FIG. 6: Energetic surfaces of reactive system H + H2 in the coordinates plane (r, 9) for different 
fixing coordinates R. It is evident from pictures in the strong interaction region (the region where 
interactions between all particles are essentially) that there are two potential holes one around the 
angle 9 = (or collinear configuration A + (BC) ) and the second around the angle 9 = ir (or 
collinear configuration A + (CB)). It is obviously that during the evolution the three body system 
with big probabilities can fall within first or second hole. 




FIG. 7: Energetic surfaces of reactive system H + H2 in the coordinates plane (R,8) for different 
fixing coordinates r. On these energetic surfaces are generated gradients which push on the three 
body system to collinear configurations. 

APPENDIX B: MASS-SCALED DISTANCES BETWEEN THREE-BODY 

For right introducing of NCC system it is necessary to connect the curve of reaction 
coordinate Qfj/ in the plan of (qo,Qi) with two equilibrium distances of diatoms of (in) and 
(out) asymptotic channels. Let us consider the distant between A and B particles in the 



(in) Jacobi coordinates (see FIG. 1): 

Rab = \/R 2 a - 2pR a r a cos 9 a + (pr Q ) 2 , p = — ^ — , (Bl) 

m b + m c 

where indexes which note channels are admitted. After simple transformations taking into 
account ([2]) and ([3]) from (IBlj) we can find: 

Rab = y^-^o - 2 P<?o<?i cos 9 + \ 2 (pqi) 2 = A" 1 ^ 2 - 2bq oqi cos 9+ (b Ql ) 2 = \~%, (B2) 

where b = p\ 2 = cot <p, and q<i is mass-scaled distance between A and B particles. Now we 
can investigate the asymptotic behavior of mass-scaled distance qi- From Eq. ( 1B2I) with 
taking into account (jSJ) in the (in) channel (in the limit of ql — > q~ ) may be find: 

Q2\^ f -q C o= (q c 2 -q c )\ qC ^ q -=-bq-cos9. (B3) 

From which is following that curve 3^ in the reactant (in) channel comes to equilibrium 
diatom (BC) distance q~ . 

Because Jacobi scattering angle which is defined in (in) channel in the (out) channel is 
limited to zero ( 6 — > 0) from (1B2|) and (jSJ) can be find: 

(l2\ q ^ oo =q c -bq c 1 = qt q , (B4) 

which is mean that in product (out) channel curve 3j/ coming to equilibrium distance q£ q 
of (AB) diatom. 

So, the definition of curve of reaction coordinate 3j/ which is defined by formulas (jSJ) 
and (El) satisfies the aforementioned asymptotic conditions. 

Completely note that definition of curve of reaction coordinate 3^ for exchange reaction 
A + (CB) — > (AC) +B we can find from (j3J) putting instead of constants b and g+ constants 
(1 — 6) and q£ correspondingly, where q^ is an equilibrium distance between pair (AC) in 
the asymptotic channel (q\ — > +oo). 

As for the dissociation reaction A + (BC) — > A + B + C it can be described by term of 
both NCC systems. Only in this case we must remember that Jacobi angle does not limit 
to zero and can have any value. 



APPENDIX C: CALCULATION OF METRIC TENSOR IN CURVILINEAR RE- 
ACTION COORDINATES 



Let us calculate the transformation between Jacobi scaled and NCC coordinates systems 
in the plan (q , q{). As well known covariant metric tensor defined by: 



7y = Yl 



k=0 



dq k dq k 

dx l dxi 



{x 1 } = (x ,^ 1 ) = (u,v), i,j,k = 0, 1. 



(CI) 



In order to do this, first we calculate the non-diagonal terms. Using dSJ) and (151) , as well as 
the fact that x° and x 1 are orthogonal, one can obtain from (IC1I) the following expression: 

' dq\ 



7oi = 7io 



— v cos (p(u) — 

du du 



sin <p + 



— u sin 0(iz) — 
du du 



COS ( 



— — sin (p[u) 
du 



dqj 
du 



COS 0{U) 



Requiring orthogonality of the NCC we set 7 01 = 7 10 = 0, giving: 

cot (f)(u), 



dq\ 



(C2) 



(C3) 



on the curve Qv,. Note that the function 4>(u) has a simple geometrical meaning: it corre- 
sponds to the angle between the tangential to the reaction coordinate and q\. From (1C3P it 
is easy to get: 

. s 1 dql/dql 



Slllffl lu] 



COS 0(U I 



After derivation sin0(u) by variable u with taking into account (1C4I) one can write: 



(C4) 



(sin (/>(«)) = cos (f)(u) 



dqljdql 



du 



d fdq§ 



1 + (dgg/dgf) ) 
Using equations (jSJ) and (jSJ) from (IC5I) may be find: 



2"! 3 / 2 (in \ dq\ 



(C5) 



[sin0(ti)) 



2a 



&~a/(gi ~g eg ) 2 / 1 



b - a/ (q<[ - q eq ) 4 



-3/2 



(C6) 



where symbol () ;s is denoted derivation. 

The expression for function (sin0(u)) correspondingly have a form: 



(COS0(W)) :M 



2a 



1 + 



b-a/(ql-q eq )' 



-3/2 



(C7) 



For matrix element 700 from ACID may be find expression: 



7oo 



dqo \ 

— — + = — v cos <p{u)— + — v sin <p{u)— . (Co) 



<9u 



9gi 



<9u 



du 



du 



dq{ 



du 



du 



Substituting the derivatives (IC6I) and flC7j) into (1C8|) we obtain: 



^[i + a>h*(|) : 



(C9) 



where .ff (it) is Gaussian curvature of reaction coordinate curve 3^/ and s is correspondingly 
the length along the 9^: 



if(u) 



2a 



- 



du J 



b-aj (q{ - q eq ) 



b-aj {q{ - q eq ) 

b + a/(q c 1 -q-) 



-3/2 



(CIO) 



Note that the Lame coefficient is given by 77 = y/^yoo- 

We also note that the elements 702 = 702 — 712 = 721 vanish. Remind that the second 
diagonal element 7n may be obtained by using (jSJ)-© an d fICll) : 



7n 



dq 
dv 



dqi 
dv 



(Cll) 



Thus the metric tensor 7^ is diagonal. Finally we note that the covariant tensor is defined 
as: 

Y J = B i:j />y, 7 = Det^u) = 7 00 7n = 700, (C12) 
where Bij is the algebraic complement to 7^. 



APPENDIX D: 5-MATRIX CONTRACTION BY REACTION COORDINATE 
(INTRINSIC TIME) FORMALISM 

Between full wavefunctions in the different body-fixed systems the following transforma- 
tion may be written (see (124)) ): 

^Pa 1/2 ^R(^v w e a )D J RK (n a ) = ^pT^r^^^Mk^)- (di) 

K K 

Next we multiply the equation fIDll) by the Wigner D^*-.(f2)-function and integrate, using 
the volume element dfl = sin£ 2 dC^id^dCz- With help of formulas: 

D M*M[(Cl: C2, C3)D^ 2M ^(Cl, C2, C3) dVt = ^J— ^ Sj 1 J 2 Sm 1 M 2 Sm[M^ 



and (125]) we get: 



^ K 

Now we are multiplying the equation ( Ti~9l) on the (out) asymptotic state (I21|) -( l22l) and 
integrating it by coordinates v and 9 after which in the limit u — > +00 the scattering matrix 
elements may be found: 

lim « aj 0«)<&S#& (u* Op) 



nj K 

= E ^ n '™ S j '1 S K'K S ijK <- njK = S n'j ' K> «- njK ■ ( D3 ) 
nj K 

So we found a new expression for the scattering S-matrix elements between (in) and 
(out) asymptotic states which are defined on different a and (3 surfaces: 

Si>j>K>^njK= km (^^(u a ,v a ,9 a )^%(up,vp,9p)) (D4) 

Note that in this approach coordinate of translational motion u plays similar role as a 
usual time in the standard scattering theory and will be called intrinsic time. 

Now it is important to find the S'-matrix elements form in the one coordinates system a 
or correspondingly /3. 

Putting the expression for overall wavefanction fl27j) in the (1D4I) we can find: 



S J nTK ^n 3 AE) = lim J^(E4,W<1 J («.M)^(«^^)) . (D5) 



«^+oo y Pnj ' n 3 ....... . /vfi0f) 

Since the angle $ between vectors q 0a , and q 0/3 in the (out) channel (when u — > +00) is 
limited to ir and taking into account that: 

d(K-K)^) = (-lf- K S KK , (D6) 
it is easy to get the expression for 5-matrix elements: 



Sl, jlK ^ njK (E) = ,,^( $ £# ( w > v, 9)$ { °pK>(u, v, 0)) a . (D7) 



Besides in the expression ( ID 71) the coefficient ^Jp^/p n j is described flux normalization con- 
stant: 

lr> + , ., / n \ 1/2 

(D8) 



lim ( — ^ 



Pnj u^+oo\ Pj3 . 

In another words we proved that it is possible to represent the full and asymptotic wave- 
functions in term of same NCC system without matching between reactant and product 
channels. It is very important result of developed theory. 

So, the matrix elements S n ij'K'^njK(E) describe the amplitudes of transition probabili- 
ties between the sets of quantum numbers (njK) and (n'j 'K') correspondingly in the (in) 
and (out) channels at fixed total energy E and will be called reactive S'-matrix elements. 

APPENDIX E: SOLUTION OF CLASSICAL OSCILLATOR PROBLEM 



Now we turn to investigation of equation ( 15TT) with natural boundary conditions ( |53l ). 
This equation can be solved with the help of hypergeometric functions. The solution of Eq. 
(I5T]) . which has the form of an incident wave with positive frequency uoi > in the (in) 
channel (when u — > — oo) for a model frequency fl52|) may be written in the form (see for 
example 54|): 



uj. 



A ' 



where 



— i- 



OirAT) 



UJ; 



■1/2 



CJ_ UJ; 1 / ,, .« 

1 + -(l + tanh(Ar)) 
A 2 



exp |zu; + t + i — In [2 cosh (Ar 

-1/2 /. \ 

— cc 0J i exp (lUJiT), 



(El) 



(E2) 



Now we can write the solution of Eq. (I5TT) which satisfies another asymptotic condition in 
the (out) channel: 



UJ 



f 



-1/2 



expjzcj+r + i— In [2 cosh (Ar)]| x 



1 - -i°^\ 1 - z'^y-; -(l - tanh (At)) r ^+oo uj f 1/2 exp (iuJ f r) . (E3) 



A ' A ' 2 

Both these solutions are connected with each other by Bogoliubov 



551 ] linear transformations: 



(Tin(r) = Cia out (r) - C 2 a* ut (r), 



(E4) 



Coefficients C\ and C 2 are easily calculated, using of properties of hypergeometric functions: 

'u^N 1 / 2 r(l +iui/X)T{iu f /X) 



Cl =< <Jin(r)(J out (T) > T = ( ) 



uji ) T(l + iuo + / X)T(iuo^/ X) 



C 2 = - < a^aUr) >r= {-) r(l - ^ + /A)r(-^_/A) " (E5) 



Recall that similar to (IE4I) . the following expression may be written for the solution <y ou t{ r )'- 



CTout(r) = C 1( r in (T) - C 2 <(r), (E6) 
Now taking into account (1E4I) . (1E5I) and (1E6I) for coefficient Ci and C 2 we obtain: 



2 = sinh 2 (7r^ + /A) 2 = sinh 2 (7r^/A) 

sinh(7ru; in ) sinh(7ru; oltf ) ' 2 sinh(7ru; iri ) sinh(7ru; out ) ' 



Finally from ( 1E7|) we can find the normalization expression: 

|d] 2 — ]C 2 | 2 = 1- (E8) 

So, taking into account ( 1E3I) and ( 1E7I) . we get the exact solution of Eq. (JED), which we must 
use for future analytical computations, in the form (IE4I) . 
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